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Abstract: Prospects for reaching persistent entanglement between two spatially separated 
atomic Bose-Einstein condensates are outlined. The system set-up comprises of two 
condensates loaded in an optical lattice, which, in return, is confined within a high-Q optical 
resonator. The system is driven by an external laser that illuminates the atoms such that 
photons can scatter into the cavity. In the superradiant phase a cavity field is established 
and we show that the emerging cavity mediated interactions between the two condensates is 
capable of entangling them despite photon losses. This macroscopic atomic entanglement is 
sustained throughout the time-evolution apart from occasions of sudden deaths/births. Using 
an auxiliary photon mode and coupling it to a collective quadrature of the two condensates we 
demonstrate that the auxiliary mode’s squeezing is proportional to the atomic entanglement 
and as such it can serve as a probe field of the macroscopic entanglement. 

Keywords: cold atoms; condensate; entanglement; cavity 


1. Introduction 

After years of experimental activity in cavity QED mainly exploring fundamentals of light-matter 
interactions at a true quantum level [1], in the second half of the 90’s, these types of experiments took 
a new direction; it was demonstrated that they may constitute a platform for implementing quantum 
information processing [2]. In [3], for example, the light-matter interaction was monitored such that the 
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photon state became entangled with a single two-level atom and then by letting a second atom ‘absorb’ 
the light field a two-atom entangled state was obtained. In this protocol, after the first step quantum 
information is stored in the light field before it gets transmitted to the second atom. Such a scheme is 
inevitably sensitive to photon losses in the middle step, which, therefor, requires short operation times. 
It was, however, realized that the ‘cavity mediated atom-atom interaction’ could also be utilized for 
directly entangling two atomic qubits without ever populating the lossy photon mode [4]. In such a 
scenario, photons are only virtually excited, but nevertheless, integrating out the photon field results in 
an effective atom-atom interaction. Then, by carefully tuning the interaction time the desired entangling 
operation could be implemented. Somewhat surprisingly, even for a ‘classical’ photon field, i.e. thermal 
state, it is indeed possible to use it for entanglement generation between the two atoms [5], and thus one 
may expect that the entanglement preparation is prone to both photon losses and other possible noise 
sources originating from the quantized radiation field. 

Following the pioneering cavity QED experiments, it has in recent years become possible to 
coherently couple not only a single atom but a gas of thousands of ultracold, condensed, bosonic atoms 
to a single quantized light mode [6] . This scales up the effective light-matter interaction strength, and 
the strong coupling regime has indeed been demonstrated in several groups^ [7]. At these sub-Kelvin 
atomic temperatures, the atomic motion must be treated quantum mechanically - it couples to the 
light field via the photon recoils. This led to explorations of a new, ultracold, regime of cavity QED. 
Much of the focus has especially been devoted to the collective phenomenon of superradiance^ [8], 
realized by coupling the light degrees-of-freedom to the atomic vibrational modes rather than to internal 
atomic electronic states [6]. Here, an atomic condensate is trapped inside the resonator and an external 
‘perpendicular’ drive/pump scatters photons into the cavity by illuminating the atoms. It turns out that at 
a critical pump amplitude, the atoms self-organize themselves in such a way that the different scattering 
processes constructively interfere making it possible for the photons to enter the resonator [9]. Below 
the critical point, the condensate is approximately uniformly distributed and scattering processes add up 
destructively leading instead to a vanishing cavity field. This change in the system properties has shown 
to share many resemblances with the Dicke phase transition from a normal (vacuum cavity field) to a 
symmetry broken superradiant phase (coherent cavity light field) [10]. Properties of this transition has 
been the subject of several experimental activities^, for example the symmetry breaking has carefully 
been mapped out and also the critical exponents [1 1,14]. 

Eike for the single/few atom cavity QED realizations, a natural question arises whether also these 
‘many-body’ cavity QED systems can serve as toolboxes for quantum information processing. A first 
step along these lines would be to study possibilities of entanglement generation among the many 
atoms. A crucial difference compared to the traditional cavity QED schemes is that these systems are 


* For N two-level atoms equally coupled to the field, the collective coupling goes as VN, and thus the characteristic time 
scale (inverse of the Rabi frequency) is greatly reduced. 

^ In this setting, superradiance does not necessarily manifest itself as an enhanced radiation process [8] but instead as 
collective properties of the atomic condensate. 

^ An alternative dynamical realization of the Dicke type transition is to consider internal electronic states Raman coupled 
via an external laser, first proposed in [12] and later experimentally implemented in Ref. [13]. 
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driven. Thus, one should rather eonsider possible persistent entanglement in the long time limits of 
the driven-lossy system [15]. Just like for Refs. [3-5], the photon mode plays the role of an aneilla 
system mediating atom-atom entanglement. The situation is different from those of [4,5] though sinee 
the non-vanishing field amplitude is sustained from balaneing the gain from the input field to the losses 
due to photon deeay, while in Refs. [4,5] the system is elosed and the photon field initialized in some 
state. Furthermore, in the present setup the atomie fields will inevitably get entangled with the photon 
field, whieh eauses an effeetive deeoherenee for the atomie state, i.e. the redueed state for the atoms will 
be mixed. In addition, photon dissipation gives rise to another deeoherenee ehannel for the atomie fields. 
All in all, even if in the few atom seenario quantum eoherenees of the radiation field is not neeessary for 
atomie entanglement generation [5], due to the differenees with the present many atomie situation it is a 
priori not elear whether any atom-atom entanglement may survive. 

Inspired by a new system setup of the Esslinger group [16] we analyze the above problem by 
eonsidering two spatially separated atomie Bose-Einstein eondensates (BECs) loeated within a eavity 
and perpendieularly driven by an external laser. Entanglement properties between the two atomie 
eondensates (as previously mentioned the atomie degrees-of-freedom are eolleetive atomie vibrational 
modes) are studied. As eaeh eondensate eonstitutes a maeroseopie meehanieal oseillator, the BEC-BEC 
entanglement is manifested in the eorresponding vibrational states. By ealeulating the logarithmie 
negativity, persistent entanglement is demonstrated even in the open system ease when photon losses 
are taken into aeeount. More preeisely, the BEC-BEC entanglement displays an infinite series of 
entanglement sudden deaths/births sueh that on average the two eondensates keep their entanglement 
for an infinite time. Entanglement between the two eondensates is, however, only possible in the 
superradiant regime. Close to and at the eritieal point, the main part of the quantum eorrelations are 
shared between the photon field and the atoms, while deeper in the superradiant phase the BEC-BEC 
entanglement beeomes more prominent. This is an example of entanglement sharing in a tripartite 
system [17]. In order to probe the entanglement between the two eondensates we introduee an auxiliary 
photon mode and eouple it to a eolleetive quadrature of the two eondensates. We demonstrate that near 
the eritieal point the squeezing of this auxiliary mode seales linearly with the the BEC-BEC entanglement 
and thus provides an indireet readout of the entanglement. However, losses of this additional photon 
mode implies that in the long time limit no quantum entanglement survives between the two eondensates. 

2. Model system 

By now, aehieving a eoherent light-matter eoupling between a single BEC and a quantized photon 
mode has been demonstrated in several groups [6,7,11,18,19]. The eoupling between single photons 
to the atomie eolleetive vibrational modes realizes a longstanding goal of optomeehanies, i.e. eoherent 
interplay between photons and a maeroseopie meehanieal oseillator prepared in its ground state [20]. 
This has spur great aetivity in exploring hybrid quantum systems with the objeetive to devise sealable 
quantum arehiteetures [21]. In partieular, generating nonelassieal states in atom-photon eoupled hybrid 
quantum systems has reeeived signifieant theoretieal and experimental interest [22,23]. Continuing this 
quest, we envision two spatially separated BECs eonfined inside an optieal resonator and explore whether 
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macroscopic entanglement between two atomie BECs can be generated via the eoupling to a eommon 
photon mode. 

With the partieular system in mind, the model Hamiltonian is given by [9] 


H = -Aah 


+ 5 : 

i=i,2 • 


dx^J(x) 




—+ 17od^a cos^(kx) + ir^cosfkx) (d^ — d) 
2m ^ ' 


$,(x) 


( 1 ) 


Here, the parameter A is the eavity-pump detuning, L is the eavity length, m the atomie mass, k 
(= (0,0, fc)) the wave number, Uq = g^/Aa is the effective atom-eavity eoupling with go the bare 
Rabi frequeney and the atom-pump detuning, and g = ilgo/Aa with H the pump amplitude. The 
photon ereation (annihilation) operator is d^ (d), and ^j(x) is the atomie field operator annihilating an 
atom in eondensate j at position x. The operators obey [d, d^'] = 1, ^i(x'),^j(x) = 6ij, and the 

remaining eommutators are identieally zero. In deriving (1) the exeited atomie eleetronie state has been 
adiabatieally eliminated under the dispersive assumption |A|, |Aa| S> go, and furthermore we consider 
a vanishing overlap between the two eondensates sueh that the eross terms of the atomie operators ean 
be negleeted. 

As a next step, we seeond quantize the Hamiltonian by expanding the field operators in the two lowest 
vibrational modes [9,1 1] 


^.(x) = + 


— cos(kx)c;^ 


,(i) 


( 2 ) 


where annihilates an atom of vibrational mode 0,1 in eondensate j. Inserting the operators (2) into 
the expression (I) for the Hamiltonian gives the low energy seeond quantized model, now eontaining 
five independent boson modes. Taking the eonservation of number of atoms of the two eondensates into 
aeeount, it is praetieal to employ the Sehwinger spin-boson mapping [24] 


Sx 2 ^ 


(i)t s(i) , 

Cq -t- Cq Cl 


T = 

-L rf. 


_ . ^(2)t-(2) , 42)t^(2) 




(l)tdl) 


Cq Cq c 


(l)tA(l) 
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. _ A(2)t^(2) 
Cq Cq Cl 


( 3 ) 


^. = 2 ("1 


(i)tAi) A(i)tAi) 




T, = -(cr'c( 


(2)]42) 42)]42) 


*^0 ^0 


to reeast the Hamiltonian in a generalized Dieke model 

H = —Ad^d -f ujr^Sz -f Tj) -|—^=(d^ -f- d) (>S-|_ -f S—) — -j= (d'^ -f d) (T_|_ -f T_) -f —a^a{Sz-\-Tz), (4) 

y N y N A 


where ur = hk‘^/2m is the reeoil frequeney. The atomic operators {S±, Sz,f±,Tz} satisfy the SU(2) 
eommutation relations [J±, Jz] = tJ±, [J+, J-] = “^Jz, where J G {5,T}. The standard Dieke model 
is laeking the last term and the boson mode eouples only to a single spin We have introdueed 
the total number of atoms N ; the effeetive atom-light eoupling g = \/2Ng and the Stark shift term 
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Since the spins are preserved we could define a collective spin R = S+T and thereby recover the regular Dieke spin-boson 
interaction. 
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U = NUq/A. With this parameter sealing, for large atom numbers every term in the Hamiltonian 
0{N). This ensures that in the thermodynamie limit, N, L ^ oo and N/L =eonstant, the model 
exhibits a seeond order phase transition [10]. In the following we will omit the Stark shift term, i.e. 
setting U = 0, valid in the regime were the pump amplitude is mueh larger than the Rabi frequeney 
whieh is justified in typieal experiments [11]. In absenee of any losses the eoupling Qc = 
separates the two phases, normal {g < g^) from the superradiant {g > gc) one. In the normal phase, 
the ground state is simply the vaeuum state of all modes {Sz = Tz = — A^/4). In the superradiant 
phase, however, a Z 2 parity symmetry is spontaneously broken whieh is for instanee refleeted in a 0 or 
TT phase of the photon field [14]. The spins in this phase are tilted away from the south pole toward 
the ‘x’-direetion [25]. Classieally, the phase transition manifest as a pitehfork bifureation [26], and 
it has been argued that entanglement generation should be most pronouneed in the symmetry broken 
phase [27]. The eavity indueed atom-atom interaetion whieh underpin the entanglement is transparent 
when the eavity field is eliminated giving rise to an effeetive atomie model of the Lipkin-Meshkov-Gliek 
type [28]. Entanglement properties of sueh models have been studied in the Refs. [22,28]. Integrating 
out the photon field and imposing an adiabatie approximation in order to derive a Lipkin-Meshkov-Gliek 
model is not always justified, for example for moderate detunings A and good eavities (low photon deeay 
rates k). Clearly, in this ease the full system has to be eonsidered. This may be extra important when 
studying entanglement properties sinee the presenee of the eavity field will itself get entangled with the 
atomie fields and thereby indireetly degrade atom-atom entanglement. In addition, the photons may leak 
out of the eavity whieh result in an additional deeoherenee meehanism for the atoms. 

To take photon deeay into aeeount we follow the standard proeedure by eoupling the eavity mode 
to a set of harmonie bath oseillators and apply the Bom-Markov and seeular approximations that are 
supposed to be valid in the optieal regime [29] . Doing so we end up with a Markovian master equation 
of the Lindblad form. Assuming a zero temperature bath the resulting open dynamies ean be modeled 
through the following master equation 

^p=-i[H,p] +nCalp], ( 5 ) 

where jCx[p] = 2xpx^ — x^xp — px^x is a Lindblad superoperator [29] and p is the full system density 
operator. We should keep in mind that we eonsider a driven system out of equilibrium. This justifies 
a phenomenologieal introduetion of losses as modeled through the master equation (5). In partieular, 
the light-matter eoupled Hamiltonian (4) is written in the frame of an external drive imparting it a 
non-equilibrium eharaeter. Lor a time-independent system eoupled to an external bath it is required 
that the steady state will obey the prineiples of equilibrium statistieal meehanies, while here no sueh 
restrietion is enforeed on the dynamies whieh arises out of an implieit time-dependent Hamiltonian [30] . 

2.1. Semiclassical analysis 

As we are interested in entanglement generation on a maeroseopie level we assume the number of 
atoms N ^ 1. This is typieally the ease in experimental realizations were N > 10^ [11,14,18,31]. 
Lor partiele numbers of this size, mean-field approaehes qualitatively eaptures many of the system 
eharaeteristies [11]. However, sueh a semielassieal mean-field approaeh is ineapable of eapturing 
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or quantifying the quantum correlations. For bosonic (spin) systems, such quantum corrections are 
conveniently explored in terms of the Holstein-Primakoff representation [32], from which we can extract 
information about the entanglement. Before analyzing the quantum features we first need to determine 
the mean-field solutions of the model which are the steady state solutions of the semiclassical equations 
of motion governed by the master equation (5). 

Let us define (a) = a, {S_) = /3,{T_) = 6,{Sz) = w^, and (T^) = Wf. The semiclassical 
equations of motion are given from the Heisenberg equations under the factorization assumption, i.e., 
((a'*' + a)(J_ — J+)) —)■ ((a'*' + a))((J_ — J+)), ((a'*' + d)Jz) —)■ ((a'*' -f a)){J^). This results in a finite 
set of coupled nonlinear equations 


a 

$ 

Wg 

6 

Wf 


— {k — iA)a — -\- (3* -\- 5 + 5*), 

-h 2i-j={a + Q*)wg, 

—iuJnS + 2i-^j={a + a*)wf, 
v N 

{a + a*){6-6*). 


y/N 


( 6 ) 

(7) 

( 8 ) 
(9) 

( 10 ) 


At this level we do not include any back-action of the reservoir apart from the dissipative photon decay. 
This is justified when the cavity mode is coupled to a heat-bath in thermal equilibrium and the associated 
Langevin noise operators thus satisfy = 0 [33]. Spin conservation for each atomic BEC at the 

mean-field level takes the form 


wl + W = 

+ = X- 


( 11 ) 

( 12 ) 


Making use of the above equations we can find the steady state solutions for a, /3, 5, Wf, Wg. We find that 
the solution predicts the occurrence of two distinct phases (normal and superradiant) which are separated 


by a critical value of the light-matter coupling = 




.Note how the critical coupling is modified 


2V2^jAi 

by the photon decay rate k [9,12,34]. The open nature of the system not only shifts the transition, but 
also alters the critical exponents [35]. 

For g < gc, the steady state solution is a fully inverted (trivial) state with a = f3 = 6 = 0, Wf = Wg = 
—N/2. This represents the normal phase. For g > g^, however, this trivial state becomes unstable and 
two new stable steady states arise with 


a = =F2i 


g^ 

K, — iA 




(3 = 6 



N 

~2 


(13) 


Wf = Wg = 


9c 

9 


2 
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This is the aforementioned pitehfork bifureation [26] . In partieular, we will show later that a eonsequenee 
of the light-matter eoupling exeeeding a eritieal value {g/Qc > 1) is the generation of bi-partite 
entanglement between the two distant BEC samples, while they remain separable in the normal phase. 


2.2. Linearized light-matter interaction 

As argued above, in the large partiele regime the quantum fluetuations are small and may be treated in 
a linearized approaeh. Starting with the normal phase we use the Holstein-Primakoff representation [32] 


s_ 

= yN-bWb, 

(14) 

s. 


(15) 


= \/ A — cf cc, 

(16) 


C C 2 5 


T 

(17) 


where b and c are the two new boson modes deseribing the bosonie exeitations around the mean-fields 
(13), to reformulate the Hamiltonian (4) in terms of three boson fields. By further making use of the 
large N limit, we linearize the resulting Hamiltonian to arrive at 

= —Ad'*'d -I- Ufi{b^b -1- c'*'c) 4- g{d'^ -f a)(6'^ -f 6 -f c'*' -1- c). (18) 

In the superradiant phase (g/gc > 1) the eavity and the atomie field modes are maeroseopieally exeited 
and the expansion of the Hamiltonian (4) is preeeded by first displaeing the operators d, b, c around their 
steady state amplitudes (13). This gives the following linearized Hamiltonian 

fTg™ = — Ad'*'d -I- Vtilj^b -\- c'*'c) -f + <? + c'*'^) -1- 0(d'*' -\- d)(b^ -\-b-\- -\- c), 


(19) 


where 


2/r 4/i(l -f /r) 


c = 


a;R(l - /i)(3 -f n) 


8fi{l -f /i) 


= 9d 


d = 


1 -f /i 

2 


Note that in the limit g —)■ g^, equals In these linearized forms, both the Hamiltonians 

are quadratie for whieh entanglement properties have been thoroughly studied [36]. For example, 
the logarithmie negativity N [37] used in the next seetion is a good measure in order to quantify the 
bi-partite entanglement. To further simplify the Hamiltonians (18) and (19) we define the following set 
of eolleetive operators 

p = (6-f c)/2-f d/\/2, 
q = (6-f c)/2 - d/\/2, 

= {b-c)/V2, 


s 
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or, inversely 

a = {p-q)/V2, 
b = {p + q)/2 + s/V2, 
c = {p + q)/2 - s/V^. 


In terms of these the Hamiltonian in the normal phase (18) takes the form 


= ujpp^p + ojqq^ + ojrs'^s + g{p^q + q'^p) + -^{p^ + p^"^ - - q^"^) 

v2 


( 20 ) 


with 


Up 


Uq = 


9 = 


Ur- A 

+ V 2 g, 

2 

Ur — A 

-V 2 g, 

2 

Ur + A 


2 



( 21 ) 


Similarly, the Hamiltonian in the superradiant phase (19) ean be expressed as 

= Upp^p + Uqq^ + {p^q + q'^p) + C{pq + q^p'^) + C(g^ + 


where 


Up 


Uq = 


.2 V2 
Q-A 
H- A 


( 22 ) 


720, 

720 . 


From the above expressions it is seen that the relative mode s deeouples from the remaining two modes. 

So far, the open eharaeter of the system only enters in the system parameters of the above effeetive 
low energy models. We have not taken the Lindblad terms into aeeount. Sinee losses only appear for the 
photon mode, the Lindblad term of Eq. (5) is unaltered going to the Holstein-Primakoff representation. 
Nevertheless, working in the eolleetive basis (20) we should also transform the loss term; 

d '' tv tv 

—p = -i[H, p] + -Cp[p] + -Cq[p] - -{2ppq^ - q^pp - pq^p) - -{2qpp^ - p^p - pp^g),(23) 


where H = or H = depending on the phase studied. Note that sinee the s mode is purely 
atomie it does not eouple to the other modes even when the system internets with the surrounding 
reservoir. In partieular, in the normal phase of the elosed system the ground state for the s mode is the 
vaeuum. In the superradiant phase, however, the eorresponding ground state is a squeezed vaeuum [8]. 
Thus, if the system’s s mode is initialized in vaeuum it is not an eigenstate and in partieular it evolves 
aeeordingly 
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A consequence of this time-dependence of the s mode is that correlators involving the modes b and c 
may also become time-dependent, even though the collective modes p and q may approach a stable fixed 
point. For instance, if the d, b, and c modes are initialized in their respective vacuum states then in the 
steady state one gets 

1 idC^sin^ 

(6^6)ss = (c^c)ss = ■^{p^p + q^ + p^ + q^p)ss - ^ _ ^2 -(^5) 

Similar expressions can be obtained for other correlators containing the b and c operators. 

3. Macroscopic entanglement 


3.1. Sustainable entanglement generation 

It should be remarked that the master equations (23) preserve the initial Gaussian character of the 
modes p, q, and s (or, equivalently of modes d, b, and c). Using standard quantum optical techniques 
it allows us to convert (23) into partial differential equations for the quantum characteristic function 
defined as xi^p^ ^q) = We will not enter into the details of this exercise, but will 

refer the reader to the same approach described in [30]. As pointed out earlier, the solutions verify that 
the b and c modes remain separable throughout the normal phase. In the superradiant phase, however, 
the atomic b and c modes can become entangled^. We characterize the bi-partite entanglement between 
the various modes in terms of the logarithmic negativity N[37]. 

For times t = — 4(^, fc G Z the b and c modes do not directly couple to the s mode. 

Bi-partite entanglement between the b and c modes for such a choice of t is shown in Fig. 1 (b) and 
(d). Also shown in Fig. 1 (a) and (c) is the bi-partite entanglement between the cavity mode and one 
of the atomic BECs. As can be seen from Fig. 1 (a) and (c), bi-partite entanglement develops between 
the cavity field and each atomic cloud. This is true in the normal and superradiant phase. However, 
it is only when g/gc. > 1 that the entanglement between distant atomic BECs develops. Entanglement 
between the cavity field and each atomic cloud is maximal near the critical point g/gc = 1, which 
may serve as an indicator of the quantum critical point. This is indeed a universal behavior; at the 
critical point where the length scales (in system with short range interactions) diverge the entanglement 
also diverges [38]. Eurthermore, as a function of the effective detuning 5 = —A — ojr the maximum 
entanglement is established at the critical point and at resonance 5 = 0. It is interesting to observe that in 
our model, which does not belong to the class of short range interaction ones, the bi-partite entanglement 
reaches its maximum at the critical points only in the cavity-BEC subsystems, but not in the BEC-BEC 
subsystems. In fact, as bi-partite BEC-BEC entanglement builds up (increasing the coupling strength) 
the cavity-BEC entanglement decreases. This is the phenomenon of entanglement sharing [17], which 
for a tripartite system implies that simultaneous large bi-partite entanglement among all constituents is 


5 


We also find that the a and b modes (or, equivalently, modes a and c) do exhibit bi-partite entanglement both in the normal 
and superradiant phases. 
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< 



Figure 1. Bi-partite entanglement (logarithmic negativity N) between the cavity mode and 
each atomic BEC, (a) and (c), and bi-partite entanglement between the two atomic BECs, (b) 
and (d). The fact that the cavity-BEC or the BEC-BEC entanglement roughly complement 
one another in terms of amplitude is an outcome of entanglement sharing [17]. We also note 
that atom-atom entanglement occurs only in the superradiant phase and beyond a critical 
value of the detuning, i.e. | A| > Ac, i.e. it is not symmetric with respect to the photon-atom 
detuning 5 = —A — ur. In (a) and (b) A = —2, while in (c) and (d) k = 0.05 and for all four 
plots t = kn /fc G Z. The value of the loss rate k has been taken small, but still 
experimentally relevant [11,14,31], in order to achieve as large entanglement as possible. 
The (scaled) recoil frequency ujr = 1. 


not allowed. We note that for a special case of our model this effect has been demonstrated namely for 
the closed Tavis-Cummings model {i.e. the Dicke model with the rotating wave approximation) for two 
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two-level atoms [39]. An additional feature evident from Fig. 1 is that BEC-BEC entanglement appears 
only beyond a eritieal value of the detuning, i.e. A < Ac. In other words, being in the superradiant 
regime is not always suffieient to warrant BEC-BEC entanglement. 

The results reported in Fig. 1 are obtained for a partieular ehoiee of parameter t. Experimentally 
it is more relevant to ask the question how the time-averaged entanglement depends on the system 
parameters. One may, for example, imagine a eontinuous measurement and then averaging the results, 
or making measurements at random times and then averaging. Thus, we evaluate the time-averaged 
logarithmie negativity as a measure of bi-partite entanglement among photon and atomie modes, i.e. 
we evaluate N = 'N{t)dt/T for large sampling times T’s. The resulting time-averaged logarithmie 
negativity is plotted in Fig. 2. As is evident from Fig. 2 (a), the time-averaged bi-partite entanglement 
between the two BEC samples is indeed very fragile to eavity damping. As illustrated in Fig. 2 (b), the 
time-averaged bi-partite entanglement between the two BECs is a nontrivial funetion of the eavity field 
detuning A and the eoupling strength g/gc. In partieular, the two BECs may beeome entangled also for 
red detunings A > —ujr. Naively, one eould expeet that 5 = —A — cu/j = 0 should generate maximum 
BEC-BEC entanglement sinee then the effeetive atom-atom eoupling is the largest. However, as elear 
from Fig. 2 (b) the pieture is more eomplex; a strong atom-field eoupling implies also strong baek-aetion 
of the field onto the atoms (and also an inereased influenee of the photon reservoir on the atom-atom 
eoherenee), whieh ean eounteraet build-up of BEC-BEC entanglement. This tradeoff between inereasing 
the effeetive atom-atom interaetion and deereasing the influenee of the eavity field and its deeoherenee is 
also refleeted in the non-monotonous behavior of the negativity as a funetion of g/gc. As shown in Fig. 2 
(e) due to entanglement sharing eonstraints in a tri-partite quantum system, the BEC-BEC entanglement 
grows at the expense of the eavity-BEC entanglement. 

3.2. Inferring entanglement 

A ehallenging aspeet of any seheme involving maeroseopie entanglement generation is the aetual 
experimental deteetion of entanglement. In this seetion we will propose a strategy to indireetly infer 
the degree of entanglement generated between the two atomic BECs. The scheme can be used to 
indirectly infer the degree of entanglement generated between the two BECs in the superradiant phase 
for small values of the critical coupling ratio g/gc > 1. More precisely, we suggest to make use of 
an additional (auxiliary) optical w mode to indirectly infer the degree of entanglement between the two 
atomic samples. The auxiliary mode couples to the relative atomic s mode with a strength fh as 

^Im-aux ^ ujpP^p + UJqq^ + -f ^(pg ^t^t) 

+ (^ + + (2 ~ ^ ^ + W^s). (26) 

A discussion on an actual physical realization of the above Hamiltonian has been left as an open question 
which might be studied in a future work. It should be pointed out that the Hamiltonian (26) has a 
beam-splitter (passive) interaction between the relative s mode and the auxiliary w mode [40]. In this 
sense, the w mode initially prepared in a classical state and evolving under the Hamiltonian (26) can 
only become non-classical as a result of a quantum state swap with the s mode. The open dynamics of 
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Figure 2. Time-averaged bi-partite entanglement (logarithmie negativity N) between the 
two eondensates as a funetion of (a) k and g/Q c for A = -1, [h) A and g/Q c for k = 0.05. 
Shown in (c) is the time-averaged eavity-BEC (dashed) and BEC-BEC (solid) entanglement 
for two different values of the parameter A and k = 0.05. The rise of atomic entanglement 
as the atom-photon entanglement decreases is again an outcome of the entanglement sharing 
mechanism. Other physical parameters include = 1. 


our hybrid tri-partite quantum system, now interacting with an auxiliary w mode, is described by the 
following master equation 


J/ = -’I"' 


lin—aux 
sr 






tv tv 

- 2 - PQ^p) - 2 -p^p- pp^q) -f 


(27) 

Here we have included photon decay of the auxiliary mode with a loss rate 7. It should be pointed out 
that a damped auxiliary w mode provides an indirect dissipation channel for the s mode. In other words, 
indirect damping introduced by the w mode can allow the s mode to settle to a time-independent steady 
state. Eig. 3 shows transient bi-partite entanglement between the atomic b and c modes in the presence 
and absence of the auxiliary w mode. The time decaying envelope of the quantum entanglement between 
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the atomic modes can be attributed to the indirect damping suffered by the s mode. More precisely, 
without losses of the s mode the entanglement does not reach a vanishing steady state - the sustainable 
quantum correlations derive from the lossless s mode. With losses of the w mode included, however, 
the BEC-BEC entanglement vanishes in the long time limit. The Eig. 3 also gives a clear signature of 
entanglement sudden death and birth [41], arising from the common coupling of the two atomic BECs to 
a single cavity mode. As shown in Eig. 3, entanglement sudden death and birth occur both in the absence 
and presence of the auxiliary w mode. Similar observation of entanglement sudden death and birth has 
been reported in a scheme proposing entanglement generation among distant optomechanical systems 
[42]. 

In Eig. 4 (a) we plot the logarithmic negativity both with (4/ ^ 0) and without ('P = 0) the auxiliary 
mode vs. small values of the critical coupling ratio q/qc > 1. The linear dependence in this parameter 
regime suggests a simple one-to-one mapping between the two cases. Indeed, we find that for values 
of the critical ratio g/Q c > 1, the time-averaged quantum entanglement between the b and c modes 
in the absence of the auxiliary w mode (N^=o) is proportional to the same quantity in the presence of 
auxiliary w mode (N^^o)- Thus, the effect of the auxiliary mode on the atom-atom entanglement is to 
just scale it. This feature is captured in Eig. 4 (b). Moreover, the single mode squeezing [8] of the 
w mode is also found to behave linearly with the bi-partite entanglement between the b and c modes, i.e. 

oc as shown in Eig. 4 (c). The connection between entanglement and squeezing has long been 
known [43], and as the later can be determined by means of homodyne measurements on the auxiliary 
w mode [8], this allows the possibility of indirectly quantifying the macroscopic entanglement between 
the two BECs. It should be possible to subject the degree of single mode squeezing to a post-processing 
procedure and one can then rescale the bi-partite entanglement (N^=o) according to the linear behavior 
shown in Eig. 4 (b) and (c). 

4. Conclusions 

In this work we have explored a possibility to generate macroscopic entanglement between two BEC 
atomic samples. We considered a physical setting where two laser driven spatially separated BECs 
are placed in a high finesse optical cavity where the photons can leak out of the cavity at a rate k. 
The two BEC samples interact indirectly through the same cavity field imposing an indirect atom-atom 
interaction which generates the entanglement. In the thermodynamic limit our physical model is known 
to exhibit a Dicke-type phase transition when the light-matter coupling exceeds a critical value. This 
phase transition separates the normal phase from a superradiant phase. Starting with a semiclassical 
analysis of the steady state of our driven-dissipative model, and adding a linearized treatment of quantum 
fluctuations using the Holstein-Primakoff representation we solve the resulting master equations in the 
two phases separately. We find that the two BEC samples remain separable in the normal phase but can 
become entangled in the superradiant phase. This generation of bi-partite entanglement between the two 
BECs can be seen as a signature of our hybrid system entering the superradiant phase where quantum 
coherence persists over macroscopic length scales. In particular, by entering the superradiant phase the 
cavity-BEC entanglement peaks at the critical point and almost no BEC-BEC entanglement exists at 
this point. Deeper in the superradiant phase the phenomenon of entanglement sharing occurs where the 
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Figure 3. Time-evolution of the logarithmie negativity N, serving as a measure of 
entanglement, between the atomie b and c modes in the presenee, 'h 7 ^ 0 , (green) and 
absence, = 0 (blue) of the auxiliary w mode. Other parameters include A = —1, 

Ur = 1, K, = 0.05, 7 = 0.05 k , and g/gc = 1.05. In the presence of the auxiliary w mode 
the logarithmic negativity N monotonically decreases with time. The inset shows transient 
oscillations of N,i,=o and featuring multiple occurrences of entanglement sudden death 
and sudden birth. At these short times, the influence of the decay of the w mode is not 
apparent; we note that the characteristic time scale for this decay 7 “^ = 400. Note that with 
7 = K = 0.05 we would still observe entanglement for f < 50 and thereby also several 
deaths/births of the entanglement. 

cavity-BEC entanglement is lowered while the BEC-BEC entanglement gets stronger. Even though the 
open-driven system does not approach a steady state in the general case, the time-averaged BEC-BEC 
entanglement is persistent. Since the two BECs are coupled to the same leaky photon mode, we also 
encounter an entanglement sudden death/birth evolution. We have also outlined a strategy to indirectly 
infer the degree of entanglement between the two atomic BECs. This scheme is no direct measure of the 
entanglement between the two condensates, but it should be a clear smoking gun of it. 
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Figure 4. Entanglement (logarithmie negativity) between the two atomie b and c modes as 
a funetion of g/gc (a) in the presenee ^ 0 ) and absence (th = 0 ) of the auxiliary w 
mode. The linear dependences in the two cases lead to proportionality shown in (b); N,i,^o 
vs. N^=o- Dots are numerical data, while the solid line a least square linear fit. In (c) the 
logarithmic negativity is displayed vs the single mode squeezing S'^^o of the auxiliary mode. 

As is clear, increasing the atom-atom entanglement implies an increased squeezing of the w 
mode, which suggests that squeezing would be an indirect measure of entanglement in the 
system. Other parameters include A = —1, ojr = 1, k = 0.05, and 7 = 0.05/?. 
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